Threshold of a Random Laser with Cold Atoms 



o 
o 

<N 



Luis S. Froufe-Perez^ , William Guerin^, Remi Carminati^ and Robin Kaiser^ 

^ Instituto de Ciencia de Materiales de Madrid, CSIC, Sor Juana Ines de la Cruz 3, Gantoblanco, Madrid 28049, Spain, 
^Institut Non Lineaire de Nice, CNRS and Universite de Nice Sophia- Antipolis, 1361 route des Lucioles, 06560 Valhonne, France, 
^Institut Langevin, ESP CI, CNRS UMR 7587, Laboratoire d'Optique Physique, 10 rue Vauquelin, 75231 Paris Cedex OS, France. 

(Dated: June 3, 2009) 

We address the problem of achieving an optical random laser with a cloud of cold atoms, in 
which gain and scattering are provided by the same atoms. The lasing threshold can be defined 
using the on-resonance optical thickness bo as a single critical parameter. We predict the threshold 
quantitatively, as well as power and frequency of the emitted light, using two different light transport 
models and the atomic polarizability of a strongly-pumped two-level atom. We find a critical bo on 
the order of 300, which is within reach of state-of-the-art cold-atom experiments. Interestingly, we 
find that random lasing can already occur in a regime of relatively low scattering. 

PACS numbers: 42.25.Dd,42.55.Zz 
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Random lasing occurs when the optical feedback due 
to multiple scattering in a gain medium is strong enough 
so that gain in the sample volume overcomes losses 
through the surface. Since its theoretical prediction by 
Letokhov ^i], great efforts have been made to experimen- 
tally demonstrate this effect in different kinds of systems 
[1, 0, 0, H, 01 , as well as to understand the fundamentals of 
random lasing [tI, Isj. The broad interest of this topic is 
driven by potential applications (see [9] and references 
therein) and by its connections to the subject of An- 
derson localization State-of-the-art random lasers 
(§1 are usually based on condensed matter systems, and 
feedback is provided by a disordered scattering medium, 
while gain is provided by an active material lying in the 
host medium or inside the scatterers. In general, scatter- 
ing and gain are related to different physical entities. 

Another system that can be considered for achieving 
random lasing is a cold atomic vapor, using magneto- 
optical tr aps, where radiation trapping [ll| as well as 
lasing flU, have been demonstrated. One advantage 
is the ability to easily model the microscopic response 
of the system components, which can be extremely valu- 
able to fully understand the physics of random lasers. 
However, in such system, the ability to combine gain and 
multiple scattering at the same time is not obvious, as 
both should be provided by the same atoms. The pur- 
pose of this Letter is to address this issue quantitatively. 
Note that even though new interes ting features appear 
when coherent feedback is involved we will consider 
only incoherent (intensity) feedback. 

Following Letokhov's theory, we consider a homoge- 
nous, disordered and active medium of size L. The ran- 
dom lasing threshold is governed by two characteristic 
lengths: the elastic scattering mean free path 4c [lll,^ 
and the linear gain length ig {£g < in the cases of ab- 
sorption or inelastic scattering). In the diffusive regime, 
defined as L i^cj the lasing threshold is reached when 
the unfolded path length, L^/4c, becomes larger than 
the gain length. More precisely, the threshold is given by 



i^eff > /3^V4c^g/3, where P is a numerical factor 
that depends on the geometry of the sample {P = I for a 
slab, /3 = 2 for a sphere), and Leff = r]L is the effective 
length of the sample, taking into account the extrapola- 
tion length [15]. Another important length scale is the 
extinction length, as measured by the forward transmis- 
sion of a beam through the sample, T = e~^/^^^. The 
extinction length is given by = ~ • 

Let us consider now a homogeneous atomic vapor, con- 
stituted by atoms of polarizability a{uj) at density p, 
submitted to a homogeneous pump field. The extinc- 
tion length and the scattering mean free path are related 
to their corresponding cross-section a by ^^^g^ = pcrex,sc5 
with crex(^) = koxlm[a{u)] and osd^^) = /cq/Ottx |a(a;)p 
(l8| . As we consider only quasi-resonant light, we use 
only the wave vector = uoq/c with uoq the atomic fre- 
quency. We also define a dimensionless atomic polariz- 
ability a such that a = a x GTr/Zcg (we omit the depen- 
dence on uj in the following). As it is an intrinsic param- 
eter of the cloud, it is convenient to use the on-resonance 
optical thickness 6o = pcro^, where ctq = Gtt/Zcq is the 
resonant scattering cross-section (without pump laser). 
Using these quantities, the threshold condition writes 



^^0 > 



(3lT 



Im(Q;)) 



(1) 



Moreover, we have L/4c = ^o|<5^P and 7^ = 1 + 
2^1 [L/4c + 2(/3 - 1)^] with ^ ^ 0.71 for L > 4c 0,0. 
Note that deeply in the diffusive regime {L isc), r] ^ 1. 

Eq. ([T]) is the first result of this Letter. It shows, 
in the diffusive regime, the existence of a threshold of 
random lasing as soon as the medium exhibits gain, i.e., 
lap — Im(a) > 0. This threshold is given by a critical on- 
resonance optical thickness, expressed as a function of the 
atomic polarizability only. Interestingly, the condition 
Im(Q;) < 0, corresponding to single-pass amplification 
(T > 1), is not a necessary condition. 

The previous result is general and does not depend on 
a particular pumping mechanism or atomic model. Let 



us now specify a gain model that will allow numerical 
evaluations of the lasing threshold and of the features 
of the emitted light. We shall use the simplest case of 
strongly-pumped two-level atoms, for which the normal- 
ized atomic polarizability at frequency uj can be written 
analytically (assuming a weak "probe" intensity) (2l| . 

1 1 + 4A^ 

^(^'^'^) = -2TTlA^T^ (2) 

+ - A + V2) - 1]^V(2A - i) 
^ ((5 + - A + il2){b + A + - n^{6 + i/2) ' 

In this expression, A = (cjp — a;o)/r is the normalized de- 
tuning between the pump frequency Up and the atomic 
transition uq of linewidth F, 5 = (cj— cjp)/r is the normal- 
ized detuning between the considered "probe" frequency 
and the pump, and Q is the Rabi frequency, normalized 
by r, associated with the pump-atom interaction. For a 
strong enough pumping power, this atomic polarizabil- 
ity allows for single-pass gain, when lm(a)<0. This gain 
mechanism is referred as "Mollow gain" (l3l.[2l| and corre- 
sponds to a three-photon transition (population inversion 
in the dressed-state basis). 

For each couple of pumping parameters {A, ^7}, the use 
of the polarizability ([2]) into the threshold condition ([T]) 
allows the calculation of the critical on-resonance optical 
thickness bo as a function of 6. Then, the minimum of 
bo and the corresponding 5 determine the optical thick- 
ness bocr that the cloud must overcome to allow lasing, 
and the frequency Jrl of the random laser at threshold. 
The result is presented in Fig. [T]for a spherical geometry 
{P = 2). The result for bocr is independent of the sign 
of A and we only show the region A > 0. The mini- 
mum optical thickness that allows lasing is found to be 
bocr ^ 200 and is obtained for a large range of param- 
eters, approximately along the line Q ^ 3 A. The op- 
timum laser-pump detuning is near the gain line of the 
transmission spectrum, i.e., 5rl ~ sign(A)VA^ + ft'^ (a 
small shift compared to the maximum gain condition is 
due to the additional constraint of combined gain and 
scattering). 

The obtained critical optical thickness is achievable 
with current technology [22], showing that random las- 
ing is possible in a system of cold atoms with Mollow 
gain. As this result has been obtained using the diffusion 
approximation, the condition L/4c = ^o|ck;p ^ 1 must 
be satisfied. This is not the case in the full range of ran- 
dom lasing parameters that we have found. For example, 
with A 1 and Q ^ 3, the critical optical thickness is 
almost minimum, bocr = 213, but L/i^c ^ 0.44. In this 
case, the threshold defined by Eq. ([T]) is at best unjusti- 
fied, at worst wrong. In order to identify in Fig. [T]the 
region in which the approach should be valid, we have 
hatched the area corresponding to L/4c < 3. Note that 
random lasing is still expected in this region, but for a 
larger on-resonance optical thickness, that would allow 
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FIG. 1: Threshold of random lasing based on Mollow gain 
(2TI I. calculated for each pair of pumping parameters A (de- 
tuning) and n (Rabi frequency) with Eqs. ()l|2p . Only the 
A > part is represented, (a) Critical optical thickness 6ocr 
to allow lasing. (b) Detuning of the random laser from 
the pump frequency. The black area corresponds to a for- 
bidden region (no gain). The hatched part corresponds to 
parameters for which the diffusion approximation is a priori 
not reliable. 

the fulfillment of the diffusive condition. The minimum 
optical thickness in the region of parameters compatible 
a priori with the diffusion approximation is 347, and is 
located in the vicinity of {A = 1,1] = 1.2}. 

This first evaluation demonstrates the need for a more 
refined transport model. In the following, we use the ap- 
proach introduced in Ref. [i^l, that is based on the radia- 
tive transfer equation (RTE). The RTE is a Boltzmann- 
type transport equation that has a larger range of 
validity with respect to the ratio L/isc than the diffusion 
equation [25*]. 

Letokhov's diffusive theory [J, and the RTE-based 
theory [23] of random lasing both rely on a modal expan- 
sion of the solution of the transport equation. In order to 
compare the predictions of both models, we focus on the 
slab geometry (/3 = 1) since the modal expansion of the 
RTE is well known in this case [19*1 (to our knowledge, no 
simple expansion is available for a sphere in the RTE ap- 
proach) . The modal approach consists in looking for solu- 
tions of the form ^^(z, u, t) = ^/^^^(u) exp(mz) exp(st), 
where ^(z, u, t) is the specific intensity {z is the distance 
from the slab surface and u denotes a propagation di- 
rection). For a given real k,, s{k,) and (j)^^s form a set of 
eigenvalues and eigenfunctions of the RTE. If one denotes 
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by so{Hi) the eigenvalue corresponding to the mode with 
the longest lifetime in the passive system, a laser instabil- 
ity appears when so{k) > in the presence of gain. The 
lasing threshold is defined by the condition so{tv) = 0. 
For isotropic scattering, this eigenvalue has an analytical 
expression valid for Hcisc < 7r/2 (lol. [23|: 



5o(/^)/c = - [£-^ - A^/tan(/^4c)] 



(3) 



where c is the energy velocity. For a slab of width L, the 
dominant mode corresponds to = tt/ Leff = tt/ (L + 2^isc) • 
In practice, this determination of is meaningful as long 
as ^ = 0.71 can be taken as a constant (independent on 
L), which is the case for L > 4c- This condition sets the 
limit of accuracy of the modal RTE approach. 

The diffusive result is recovered from the RTE ap- 
proach in the limit k,£sc <^ 1 [25]. A first order expan- 
sion of Eq. (|3]) yields s'"^^\k,)/c = - Ai:^4c/3, where 
the superscript (DA) stands for diffusion approximation. 
The condition SQ^^\tz = tt/Lq^) = leads to Letokhov's 
threshold, with /3 = 1. 

The comparison between the RTE and diffusive ap- 
proaches deserves two comments. Firstly, the gain con- 
tribution to so(/^) (first term in Eq. (|3])) is the same in 
both models. Secondly, the scattering contribution (sec- 
ond term in Eq. (|3])) is larger in the RTE model by a 
factor of at most 1.13 (when L ~ 4c)- Thus, the cor- 
rection introduced by the RTE model, compared to the 
diffusion approximation, is relatively small, as it corre- 
sponds to an increase of r^^ocr of at most a few percent s. 
This means that the diffusive model gives accurate results 
down to L ~ 4c, and that in cold- atom systems, random 
lasing can occur even in a regime of low scattering. 

In Fig. [21 we compare the minimum optical thick- 
ness obtained with both models for the slab geometry 
and with the diffusive model for the sphere geometry. 
To put forward the domain of validity in each case, we 
plot the results as a function of I//4c- As expected, in 
the range L > 4c the threshold predicted by the RTE 
for the slab geometry is only slightly larger than the one 
given by the diffusion approximation, so that the two 
curves can hardly be distinguished. For the sphere ge- 
ometry, we dashed the part corresponding to the domain 
where the diffusive model is a priori not reliable, i.e., 
L/4c < 3. Nevertheless, by generalizing the conclusion 
obtained with the slab geometry, we reasonably expect 
the threshold to be located between 250 and 300. 

Let us now turn to a first characterization of such 
a random laser. An important quantity to be investi- 
gated is the emitted power as a function of the pump- 
ing power. In the stationary regime (continuous pump- 
ing) we numerically solve the optical Bloch equations 
for a strongly-pumped two-level atom (without using the 
weak "probe" approximation that leads to Eq. ([2])) to 
obtain the polarizability at the lasing frequency, includ- 
ing the gain saturation induced by the random laser 
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FIG. 2: Critical optical thickness for different geometries and 
transport models. The lower, red curve corresponds to the 
slab geometry (width L), with RTE model (continuous line) 
and diffusive model (open circle). The upper, blue curve cor- 
responds to the sphere geometry (diameter L), with the diffu- 
sive model. The part with L/isc < 3 is dashed, as the model 
may not be reliable. 



intensity. Above threshold, the laser intensity in the 
medium ex |1^rlP is determined by the condition 
so(/^, I^^RLp) = {so would be positive without gain 
saturation). The obtained intensity is analogous to the 
intra-cavity intensity of a standard laser, and thus does 
not correspond to the emitted power P^^^^ - At equilib- 
rium, gain compensates losses, and P^^^'' is equal to the 
generated power, related to the gain cross-section (jg, i.e.. 
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c^gl^RLp with cTg = (70 -lm{a)). 



p(out) 

In order to know if the laser signal can be extracted 
from the background fiuorescence, it is particularly rele- 
vant to compare the emitted laser power with the pump- 
induced fiuorescence Ppiuo cx ao\Q\'^/{l ^ 4A^ + 2|l]p). 
From this, we compute the ratio 



c,(out) 
RL 



\n 



RL 



(4) 

We plot the result in Fig. [3] as a function of for 
a pump detuning A = 1. To obtain Eq. (|4j), we as- 
sume that both pump and laser intensities are homoge- 
neously distributed across the whole system. We also 
consider only the optimum random laser frequency, thus 
neglecting the spectral width of the random laser or any 
interaction between different random laser frequencies. 
Hence we neglect several effects as mode competition [i^ 
and inelastic scattering of the laser light. Nevertheless 
we think that the order or magnitude of the ratio laser- 
to-fiuorescence powers can be realistic for actual experi- 
ments, at least as long as only one mode of the laser is 
active [26]. For the chosen set of parameters, this ratio is 
more than 5% and hence laser emission should be measur- 
able. Its distinction from the pump-induced fiuorescence 
can be made by looking at the spectrum of the emitted 
light. Another interesting prediction of this model is that 
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of freedom they offer, together with the first estimates 
presented here, make us confident that the experimental 
reahzation of a cold-atom random laser is possible with 
current technology. 
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FIG. 3: Continuous line: Emitted random laser power nor- 
malized to the pump fluorescence power, as a function of the 
pump intensity. Dashed line: Normalized laser detuning 
The random medium is a spherical cloud of two level atoms 
with an on-resonance optical-thickness bo — 650. 



the laser emission frequency shifts as the pump intensity 
is increased [Fig. [3]. This corresponds to the shift of the 
maximum gain of the Mollow polarizability. 

In summary, we have established the possibility of 
achieving random lasing with cold atoms. The random 
laser threshold is described by a single critical parameter, 
the on-resonance optical thickness h^. In the particular 
case of a gain mechanism based on a strongly-pumped 
two-level atom (Mollow gain), our model predicts a crit- 
ical 6o ~ 300. Such an optical thickness is achievable 
in current cold- atoms experiments, e.g. by using crossed 
dipole traps [22]. We have also determined the basic fea- 
tures of the emitted light above threshold, showing that 
the random laser emission should be measurable. 

Another interesting result is that, due to the large gain, 
lasing can be obtained with a low feedback (low amount 
of scattering, i.e., L ~ 4c)- This regime is similar to that 
encountered in certain semiconductor lasers with a very 
poor cavity, and is different from the working regime of 
random lasers realized to date. This new regime could be 
numerically investigated by RTE-based simulations [25i] . 

Finally, let us stress that the model developed here has 
several limitations, so that the numbers should be consid- 
ered as first-order estimates. Firstly, we have considered 
monochromatic pumping, thus neglecting inelastic scat- 
tering from the pump. The inelastically-scattered pho- 
tons may have a non-ne^gible infiuence on the atomic 
response, as shown in [27|. Secondly, the RTE model 
needs to be extended to a sphere geometry, and to a 
medium with inhomogeneous density and/or pumping. 
This would require a full numerical solution of coupled 
RTEs for the pump and probe beams [23]. We also out- 
line that the case of the Mollow gain was chosen for the 
sake of simplicity, whereas other gain mechanisms might 
be more adapted for the search of random lasing, such as 
Raman gain or parametric gain [13]. Each gain mecha- 
nism has its advantages and drawbacks, but the degrees 
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